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ABSTRACT 

The Lya forest at z > 5.5 shows strong scatter in the mean transmission even when smoothed over 
very large spatial scales, > 50 Mpc//i. This has been interpreted as a signature of strongly fluctuating 
radiation fields, or patchy reionization. To test this claim, we calculate the scatter arising solely from 
density fluctuations, with a uniform ionizing background, using analytic arguments and simulations. 
This scatter alone is comparable to that observed. It rises steeply with redshift and is of order unity 
by z ~ 6, even on ~ 50 Mpc//i scales. This arises because: i) at z ~ 6, transmission spectra, which 
are sensitive mainly to rare voids, are highly biased (with a linear bias factor 6 > 4 — 5) tracers of 
underlying density fluctuations, and ii) projected power from small-scale transverse modes is aliased 
to long wavelength line-of-sight modes. Inferring patchy reionization from quasar spectra is therefore 
subtle and requires much more detailed modeling. Similarly, we expect order unity transmission 
fluctuations in the z ~ 3 Hell Lya forest from density fluctuations alone, on the scales over which 
these measurements are typically made. 

Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe; quasars 
- absorption lines 


1. INTRODUCTION 

At redshifts close to z ~ 3, the structure in the Lya 
forest has been shown to arise naturally from density 
fluctuations in the cosmic web (e.g. Miralda-Escude et 
al. 1996). At sufficiently high redshift, however, the 
structure in the Lya forest may instead largely reflect 
the topology of reionization, and/or a strongly fluctuat¬ 
ing radiation field. In high-redshift quasar spectra with 
extended opaque regions, significant gaps of substantial 
transmission occur (e.g. Becker et al. 2001, White et 
al. 2003, White et al. 2005). This has previously been 
attributed to a strongly fluctuating UV background, as 
expected at the tail end of reionization (Wyithe & Loeb 
2005, Fan et al 2005). 

Could these transmission gaps simply arise from un- 
derdense regions where the neutral hydrogen fraction is 
lower? The transmission in the z ~ 6 quasar spectra dif¬ 
fers significantly from sightline to sightline, even when 
one averages over co-moving length scales of ~ 50 — 100 
Mpc//i. Since the density variance is small over such 
large scales, one might naively expect that the reioniza¬ 
tion of the IGM must be incomplete near z ~ 6. 

In this Letter, we critically examine this naive intu¬ 
ition. Is rapidly increasing scatter in sightline to sightline 
flux transmission a good diagnostic for patchy reioniza¬ 
tion (Fan et al. 2002, Lidz et al. 2002, Sokasian et al. 
2003, Paschos & Norman 2005)? We find that in fact 
the fractional scatter in the mean transmissivity of the 
IGM will be large at high redshift, even for a completely 
uniform ionizing background. An analogous calculation 
applies to the case of the Hell Lya forest near z ~ 3. 

2. THE FLUX POWER SPECTRUM 

We adopt the usual ‘gravitational instability’ model 
of the Lya forest valid at z ~ 3 (e.g., Hui et al 1997). 


In particular, we assume an isothermal gas and a uni¬ 
form radiation field, to see if these assumptions demon¬ 
strably break down at z ~ 6. Assuming photoioniza¬ 
tion equilibrium, the Lya optical depth is t = AA^, 
where A oc (1 -|- z)'^ ®T“°-^/r, and the transmitted flux 
is F = 6“'’’. We study the fluctuations in transmitted 
flux, Sf = (F — (F))/(F), and the line of sight power 
spectrum of these fluctuations, Ppik). We will also refer 
to the effective optical depth, Teff = —ln(F). 

It is important to note that Ppik) is related to the un¬ 
derlying power spectrum of density fluctuations Psik) in 
a rather complicated way. First, gas pressure smooths 
the baryon-density distribution on small scales with re¬ 
spect to the dark matter. Second, a non-linear trans¬ 
formation maps density into optical depth (r oc A^). 
Third, one maps from real space to skewers in redshift 
space, projecting from 3D to ID and incorporating the 
effect of peculiar velocities. Fourth, one convolves the 
optical depth distribution with the thermal broadening 
kernel. Finally, a second non-linear transformation maps 
T into transmitted flux (F = e”’’’). These issues are well 
studied in the context of the z ~ 3 Lya forest (e.g. Croft 
et al. 2002) but are under-appreciated in the reioniza¬ 
tion literature. Previous neglect of this physics has led 
to erroneous conclusions. 


2.1. Analytic Estimates 

We will therefore require numerical simulations to 
model the flux power spectrum in detail, and its de¬ 
pendence on redshift and the mean transmissivity of the 
IGM. We can nonetheless anticipate the results with an¬ 
alytic arguments. Specifically, we estimate the point- 
to-point flux variance as a function of A and Teff. We 
perform this calculation using the Miralda-Escude et al. 
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(2001) fitting formula for the gas density PDF: 


P{A)dA = GA-^exp 


(A-2/3 - C)^ 
8^279 


dA. 


( 1 ) 


Here the parameters G and G are fixed by requiring the 
PDF to normalize to unity and to satisfy (A) = 1. The 
other parameters are given by 60 = 7.61/(1 + z) and 
b = 2.5 at z = 6. We then calculate the first two moments 
of the (unsmoothed) flux distribution, (F) = {e~^^ ), 
and (F^) = ), by integrating over the gas den¬ 

sity PDF. These estimates ignore the effect of peculiar 
velocities and thermal broadening, but we include these 
effects subsequently using our simulations. 

These integrals can be approximated using the method 
of steepest descents (Songaila & Cowie 2002), to 
give (F) ^ cA^/'^exp[—dA^-‘^], and (F^)/(F)^ ~ 

cH“3/4exp[dH°'^], where c, d, c, and d are constants, 
and A = Aj^h. We find that c = 5.3, d = 5.1, c = 0.22, 
and d = 3.5 provide a good approximation to the results 
of full numerical integrations. These scalings immedi¬ 
ately suggest that the fractional point-to-point disper¬ 
sion will be quite large when A and Teff are also large, 
even assuming a homogeneous radiation field. For in¬ 
stance, when Teff ~ 5, we expect the fractional point- 
to-point dispersion to be ap/iF) ~ 4. On the other 
hand, when Teff ~ 1, the point-to-point dispersion is only 
ap/iF) ^ 0.8. When the transmission becomes small, it 
is dominated by rare voids, which drastically increases 
the point-to-point variance. 


2.2. Numerical Calculations 

In order to make a more accurate calculation, and 
to consider transmission fluctuations smoothed over 
large scales, we calculate the flux power spectrum 
numerically. We use a Hydro-Particle-Mesh (HPM) 
(Gnedin & Hui 1998) simulation with 2 x 512^ parti¬ 
cles and 5123 niesh-points in a 40/i“3]\/[pc box, assum¬ 
ing a ACDM cosmology with (cts, d, 11^, Ha) = 
(0.84, 0.7, 0.02, 0.3, 0.7). The simulation is an implemen¬ 
tation of HPM into the parallel N-body code, Mc^ (see 
Heitmann et al. 2005 for detailed tests of Mc^, and the 
Appendix of Lidz et al. 2005, Habib et al. 2006, in prep., 
for the details and convergence studies of our implemen¬ 
tation of HPM). 

We take the baryon density and peculiar velocity 
fields from the simulation and extract artificial spec¬ 
tra in the usual manner (for more details see e.g., 
Hui, Gnedin & Zhang 1997). Specifically, we consider 
z = 4.9,5.3,5.7,6.0, and adjust A to match (F) = 
0.20,0.13,0.06, and 0.01. These values are close to, but 
slightly larger than, the observational measurements of 
Becker et al. (2001), White et al. (2003), and Fan et 
al. (2005). For each redshift, we extract 5,000 artificial 
spectra from the simulation box. 

In order to consider large-scale transmission fluctua¬ 
tions, we need to extrapolate to scales beyond the size 
of the simulation box. We do so assuming that on large 
scales, the flux power spectrum and the ID linear matter 
power spectrum differ only in normalization and not in 
shape (linear biasing, e.g., Scherrer & Weinberg 1998, 
Lidz et al. 2002), although we will show in !J31 that 
the bias increases sharply with increasing redshift. We 



k [s/km] 


Fig. 1. — Flux power spectrum as a function of scale at vari¬ 
ous redshifts, with A adjusted to the mean flux level {F) shown. 
The black dashed lines show the large-scale extrapolations of the 
simulation results described in the text. The green and red dotted 
lines correspond to the length scales L = 35, 75 Mpc//i at 2 = 6 re¬ 
spectively. Notice that fcPp(fc) decreases only slowly towards large 
scales, a consequence of aliasing, as we describe in the text. At 
2 ~ 6, the fluctuations are of order unity even on scales close to 
~ 50 Mpc/h. 


normalize the extrapolation to the simulated flux power 
spectrum at fc ~ kf /2, where kf is the wavenumber asso¬ 
ciated with the fundamental mode of the box. We gener¬ 
ally consider smoothing scales only somewhat larger than 
the size of our simulation box, Fbox = 40 Mpc/h. Fur¬ 
thermore, the flux power spectrum is relatively flat on 
these scales, and so our results are not terribly sensitive 
to the details of this extrapolation. ^ 

In Fig^ we show the flux power spectrum from our 
simulations. The amplitude of transmission fluctuations 
clearly increases rapidly with redshift, despite declining 
density fluctuations. In fact, fractional fluctuations are 
of order unity on scales of F ~ 35 Mpc/h at z ~ 6, even 
for a uniform radiation field. This implies that it may be 
difficult to discern patchy reionization from transmission 
fluctuations. We study this further in P 
Another interesting feature of Fig Q] is that the flux 
power spectrum is quite flat: kPp(k) declines only grad¬ 
ually towards large scales. This is because absorption 
spectra provide only a ID skewer through the IGM, with 
a power spectrum (ignoring redshift space distortions for 
illustrative purposes): 

r°° dk 

PUh)= ^kP,,{k). ( 2 ) 

^ For instance, if we extrapolate our results usii^a power law, 
our estimate of the transmission dispersion, Eq. O, changes by 
-- 20%. 
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Thus, the ID power spectrum falls off much more slowly 
towards large scales than the 3D power spectrum. This 
is the well-known process of aliasing found in pencil- 
beam galaxy surveys (Kaiser & Peacock 1991), whereby 
short transverse wavelength modes in projection alias 
power to long wavelength, line of sight modes. To il¬ 
lustrate this quantitatively, we consider the amplitude 
of 3D fluctuations in the non-linear baryon-density field, 
Al{k) = fc^Pb(A)/(27r^), approximating the relation be¬ 
tween the baryonic and dark matter power spectra as 
Pbik) = assuming kf = 20h Mpc“^ 

(Gnedin & Hui 1998), and using fitting formulae from 
Peacock & Dodds (1996) and Eisenstein & Hu (1999). 
We contrast this with the variance of ID baryonic fluc¬ 
tuations, A^q^,(A:||) = fc||PiD,b(fc||)/7J‘, comparing line of 
sight fluctuations of wavenumber k\\ to 3D fluctuations 

with wavenumber |fc| = A:||. The ID power spectrum is a 
factor of ~ 5.5—10.6 times larger than the 3D power spec¬ 
trum for modes with wavelength L ~ 50 — 100 Mpc/h. 
Fluctuations averaged over a skewer of a given length 
are much larger than fluctuations averaged over an en¬ 
tire sphere of comparable size. 


3. THE SIGHTLINE TO SIGHTLINE SCATTER 

Given Pp{k), the fractional sightline-to-sightline vari¬ 
ance in the transmission is (Lidz et al. 2002): 


(P)2 



sm{kL/2) 

(fcP/2) 


( 3 ) 


Shot-noise is negligible for significant stretches of spec¬ 
trum with at least moderate signal to noise (see Lidz 
et al. 2002 for details). We consider the fractional 
dispersion in the mean transmitted flux, cr^^)/(P), as 
opposed to the absolute dispersion, Note that 

Ureft = ^{f)/{F) for small fluctuations. 

We show this quantity in Fig. |21 As in the flux power 
spectrum calculation (Fig. QJ, the fractional dispersion 
increases rapidly with decreasing transmission, and falls 
off only slowly as the smoothing scale becomes large^. 
By z ~ 6, assuming (P) = 0.01 (reff = 4.6), a complete 
absorption trough (P) = 0 of length ~ 50 Mpc/h is only 
a ^ 1 — cr fluctuation, even though density fluctuations 
are small, <JppD ~ 0.2 on these scales. Transmission 
fluctuations are a highly biased tracer of underlying den¬ 
sity fluctuations at high redshift, with a bias factor of 
b = {cr(^p)/{F))/appD ~ 4 when Teg ^4 — 5 (and in¬ 
creasing at higher Teg). Although on large scales the 
flux and density power spectra have the same shape, one 
can not then simply assume that U(^p-^/{F) ~ o-ppo, as 
has previously been done (Wyithe & Loeb 2005, Fan et 
al. 2005). 

Note that state of the art reionization simulations have 
thus far been done with boxsizes of less than 10 Mpc/h 
(e.g. Gnedin 2000, Sokasian et al. 2003, Paschos & Nor¬ 
man 2005), although steps towards large scale simula¬ 
tions are being made (Kohler et al. 2005, Iliev et al. 
2005). It is common practice to generate mock spectra 
from these simulations by wrapping long lines of sight 
around these small periodic boxes many times. From Fig. 


^ Since, on large scales, the flat Pp(k) is close to a white noise 
power spectrum, oc L~^. 
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Eig. 2.— Sightline-to-sightline scatter in the mean transmission 
as a function of scale and redshift. The cyan, blue, green, and red 
curves indicate the expected scatter in the mean transmission as 
a function of scale at redshifts z = 4.9, 5.3, 5.7 and 6 respectively. 
The fractional scatter increases rapidly with decreasing transmis¬ 
sion, and dies off only slowly towards large scales. 


QJand Eq. |3 it is clear that estimates of the transmissiv¬ 
ity scatter from mock spectra generated in this way are 
underestimates. For instance, the best case scenario for 
extrapolating Paschos & Norman (2005) to larger scales 
would have a flat line of sight power spectrum on scales 
larger than their box (Pbox = 6.8 Mpc//i) and match our 
results on smaller scales. Even in this best case, we calcu¬ 
late that they underestimate the transmissivity variance 
by a factor of ~ 2 for a smoothing scale of P = 50 Mpc/h. 

Equation © does not capture the full story since the 
flux probability distribution is non-Gaussian. This is ob¬ 
vious since fluctuations in the mean transmission are of 
order unity, yet the transmission is bounded between zero 
and one. To reliably calculate confidence intervals, we 
thus need to consider the full shape of the PDF of the 
transmission, smoothed on large scales. Because the ex¬ 
trapolation of the full PDF to large scales is non-trivial 
(unlike for the simple power spectrum), we will calculate 
the PDF directly from the simulation box. Thus we ef¬ 
fectively smooth on scales of L = Pbox = 40 Mpc/h but 
underestimate the width of the PDF because of large- 
scale power missing from our box. We show the PDF of 
Teff in the bottom panel of Fig. |21 It becomes increas¬ 
ingly broad and skewed, with a long tail towards high 
optical depth at high redshift. This asymmetry implies 
a bias in measurements of Teg from samples with a small 
number of sightlines. 

The top panel of Fig. 0 shows error bands for the 
expected sightline to sightline scatter in the effective op¬ 
tical depth as a function of redshift, adopting a smooth¬ 
ing length of P = 40 Mpc/h. Even for large smoothing 













4 



4.5 5 5.5 6 6.5 


z 



Fig. 3.— Evolution of Tefj and the expected sightline to sight¬ 
line scatter. The top panel shows 68% and 95% confidence re¬ 
gions for the effective optical depth, given a smoothing length of 
L = 40 Mpc/h. Notice that the shaded regions become quite 
broad at large Teff. We assume central Tefj values of (F) = 
0.32,0.20,0.13,0.06,0.01, and 0.001, at ^ = 4.5,4.9, 5.3, 5.7, 6.0, 
and 6.5 respectively. The upper limits on Tefj at 2 = 6.5 extend 
beyond the edge of the plot. They are Tefj = 10.8,14 at 68% and 
95% confidence, respectively. The bottom panel shows examples 
of the corresponding PDFs. 


scales, the scatter increases rapidly with redshift. We 
therefore caution against over-interpreting the increas¬ 
ing scatter seen in plots of Tes{z) seen in the literature 

How does the level of scatter we predict compare with 
existing observations? We compare with recent measure¬ 
ments from Fan et al. (2005) (their Table 5), which 
indicate a dispersion of r = (2.1 ± 0.3, 2.5 ± 0.5, 2.6 ± 
0.6,3.2±0.8,4.0±0.8, 7.1±2.1) at z = (4.90-5.15,5.15- 
5.35,5.35-5.55,5.55-5.75,5.75-5.95,5.95-6.25). The 
values of the dispersion in the last two-redshift bins are 
lower limits, since some lines of sight in these bins show 
complete absorption troughs. Interpolating from our re¬ 
sults, which assume a different redshift evolution for Tgff, 

wp finrl T 19 1+0-38 +0.79 n r+0.46 +0.94 p+0.48 +0.98 

we unu T — l,^.i_0.28 -0.471 ^'0-0.32 -0.58! ^-O-o.aS -0.6l! 
o p+0.67 +1.4 t q+ 0.99 +2.0 y i+3.9 +7.1\ -ittUpt-p f+p 

'J + - 0.39 - 0 . 78 ! +-O_o .48 - 0 . 98 ! '- 1 - 0.60 -l.ob Wllcie Uie 
intervals indicate 68% and 95% confidence regions. Note 
that Fan et al. (2005) considered a smoothing scale of 
Az = 0.15, corresponding to L = 44 (56) Mpc/h at 
z = 6 (5). Our numbers assume a slightly different scale, 
L = 40 Mpc/h, throughout, due to the limited size of 
our simulation box. From this comparison, we conclude 
that our results are broadly consistent with the measure¬ 
ments of Fan et al. (2005), although we remark that the 
dispersion is not the ideal statistic to compare with since 
the PDF of Teff is asymmetric. We reiterate that our cal¬ 
culations assume a uniform radiation field and indicate a 
scatter that is already comparable to the measured values. 

Furthermore, there are several systematic effects, each 


of which likely tends to increase the observed scatter 
in the mean transmissivity relative to our predictions. 
First, as remarked previously, our simulations have a lim¬ 
ited boxsize (L = 40 Mpc/h) and therefore miss large 
scale power. Second, at high redshift the quasar con¬ 
tinuum is generally estimated using a single power-law, 
extrapolating from redward of the Lya forest (e.g. Fan 
et al. 2005). In reality, the quasar continuum has struc¬ 
ture and varies significantly from quasar to quasar (e.g. 
Bernard! et al. 2003). As a result, some of the apparent 
transmissivity fluctuations in the data may actually be 
due to fluctuations in the quasar continuum. Next, Ly¬ 
man limit systems are under-produced in simulations, yet 
they add large scale power to the transmission (e.g. Mc¬ 
Donald et al. 2005), which may lead us to underestimate 
the transmissivity scatter. Finally, metal-line absorbers 
are inevitably present in the observed Lya forest yet are 
not included in our calculations. We refer the reader to 
Tytler et al. (2004) for further quantitative discussion of 
these systematics at z ^ 2 — 3. 

3.1. Sightline-to-Sightline Scatter in the Lyf3 forest 

How about the transmissivity scatter in the Ly/3 forest? 
We can estimate the fluctuations in the z ~ 6 transmis¬ 
sivity of the Ly/3 region of a quasar spectrum using Eq. 
13 but we must consider fluctuations in the foreground 
Lya forest as well. Dijkstra et al. (2004) show that the 
total power spectrum PF,tot. is: 

Pf, tot. (k) =PF,a{k) + PF.pik) (4) 

/ dk' 

^PFAk-k')PFAk')- 

This equation reflects the fact that foreground gas, ab¬ 
sorbing in Lya, is widely separated in physical space, 
and hence independent of, gas absorbing in Ly/3. The 
last term in Eq. 0 occurs because the total correlation 
function then involves a product, which transforms into 
a convolution term. Although higher-order in the trans¬ 
missivity fluctuations (Dijkstra et al. 2004), this term 
cannot be ignored here because the fluctuations are so 
large (it increases the total variance by ~ 10-20%). 

Nonetheless, over an equivalent comoving patch, trans¬ 
mission fluctuations in the Ly/3 forest are substantially 
smaller than in the Lya forest. We can see this heuris- 
tically from iI2.1 1 {F^)/{F)‘^ ^ cA“^/^exp[(iA° '‘], and 
= Aq,/ 6.24. In practice, this works out to fractional 
fluctuations which are ~ 2 — 4 times smaller in the Ly/3 
forest. Erom the simulations, over a 50Mpc/i“^ interval, 
if we assume (Fa) = (0.2, 0.01), at z = (4.9, 6) (and thus 
{F/s)[z = 6] = 0.19 - see Oh & Furlanetto 2005 for details 
on the relation between the Lya and Ly/3 absorption in 
a clumpy IGM), we obtain '7^Ftot>/(-^tot) 0.35, com¬ 

pared to ai^F)/{F) = 0-8 for Lya fluctuations. If instead 
{Fa) = 0.001 at z = 6 (and thus {F^z = 6] = 0.085), we 
obtain (r{Ftat)!Atat) ~ 0.45, compared to (J(f)/{F) ~ 2 
for Lya fluctuations. Thus, in principle the Ly/3 forest 
is a much more robust indicator of a fluctuating radia¬ 
tion field. Unfortunately, due to contamination from the 
Lyy forest, at present the Ly/3 forest suffers from small- 
number statistics. For instance. Fan et al. (2005) have 
97 vs. 19 measurements of the Lya and Ly/3 forests re¬ 
spectively, of size Az = 0.15, spanning z = 4.8 — 6.3 in 
Lya, and z = 5.3 — 6.3 in Ly/3. The observed scatter in 
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the Ly/3 forest (their Fig. 3) is consistent with density 
fluctuations alone. 

3.2. Sightline-to-Sightline Scatter in the Hell Forest 

Finally, we point out that a very similar calculation 
applies to the sightline-to-sightline scatter in the trans¬ 
missivity of the Hell Lya forest near z ~ 3. Specifically, 
we measure the flux power spectrum in the Hell Lya 
forest using a lower redshift simulation output, adopt¬ 
ing Teff^eii ~ 4 at z = 2.9 based on the measurements 
of Smette et al. (2002). These measurements are per¬ 
formed in narrow bins with co-moving length scales of 
L ~ 18—36 Mpc/h. At z = 2.9, we find that the rms frac¬ 
tional fluctuation in the Hell Lya transmissivity on these 
smoothing scales is Heii)/(A'Heii) = 1-6 — 1.2. Again, 
we caution against interpreting the observed scatter as a 
signature of fluctuations in the Hell ionizing background. 

4. CONCLUSIONS 

We do not mean to imply that the radiation field is 
indeed uniform at z 6. Our intent is only to show 
that ruling out the null hypothesis that the scatter in 
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the high redshift Lya forest results solely from density 
fluctuations is subtle. Indeed, we find that transmissiv¬ 
ity fluctuations should be of order unity at z ~ 6 on 
scales of L ~ 50 Mpc//i, from density fluctuations alone. 
We therefore caution against over-interpreting the large 
scatter in Teg seen in the spectra of z ~ 6 quasars or the 
analogous scatter seen in the z ^ 3 Hell Lya forest. In 
future work, we intend to model the effect of inhomo¬ 
geneous reionization on the statistics of the Lya forest, 
following up on the work of Furlanetto & Oh (2005), and 
to design statistical measures to discern the presence or 
absence of these fluctuations. 
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